Whole-genome sequencing and phylogenetic analysis of coxsackievirus-A16 strains causing hand, foot and mouth disease (HFMD) in India.

Hand, foot and mouth disease (HFMD) is a common childhood infectious disease, caused by enteroviruses (EVs), which can present with typical or atypical lesions. The illness is self-limiting, but it can also have serious complications. Since 1997, HFMD infections have become endemic and have increased to epidemic proportions across the Asia Pacific region, including India. Coxsackievirus-A16 (CV-A16) outbreaks occurred in India from 2005 onwards, although the clinical symptoms were noticeably different during this period. Understanding the population dynamics of enteroviruses that cause HFMD is crucial in the post-polio era because one of the circulating strain may replace another as the dominant strain. The aim of this study is to describe the genetic features of the CV-A16 strains isolated from hand, foot and mouth disease (HFMD) patients in India. Reverse transcription PCR (RT-PCR) and cell-culture-based isolation of CV-A16 was done from the 55 clinical samples. The entire genome of the CV-A16 isolate was performed from the seven isolates. After the sequences were analysed, a phylogenetic tree was created using bioinformatics tools. The total genomic length obtained was 7411 base pairs (bp). Nucleotide similarity across various regions, including 5'UTR, P1, P2 and 3'UTR, ranged from 87.0-97.9 %, 77.0-95.4 %, 80.3-96.9 %, and 77.9-96.2 %, respectively. Correspondingly, similarities in the VP1 region's nucleotide and amino acid sequences were 91.4-96.4 % and 99.3-99.7 %, respectively. Phylogenetic analysis highlighted that CV-A16 strains identified in Pune, Maharashtra, were grouped within the same cluster. The analysed CV-A16 isolates in this study aligned with subgenotype B1c. These findings have far-reaching implications for the surveillance, prevention and management of HFMD and CV-A16. Monitoring the dynamics of CV-A16 strains, informed by the genetic characteristics identified here, will significantly impact strategies aimed at tackling HFMD and its associated public health challenges.

First detection of Varicella Zoster Virus clade 9 cases in India during mpox surveillance.

BACKGROUND: The multi-country mpox outbreak across the globe has led to the systematic surveillance of mpox cases in India. During the surveillance of mpox, we encountered cases of Varicella Zoster Virus (VZV) in suspected mpox cases amongst children & adults. This study focused on the genomic characterization of VZV in India. METHODS: A total of 331 mpox suspected cases were tested for VZV through real-time PCR, and the positive samples were subjected to next-generation sequencing to retrieve the whole genome of VZV using CLC genomics software. Phylogenetic analysis has been done in MEGA 11.0 software to identify circulating clades. RESULT: Of the 331 suspected cases, 28 cases with vesicular rashes were found to be positive for VZV. The maximum genome could be retrieved from the clinical specimens of 16 cases with coverage greater than 98% when mapped with reference strain Dumas (NC 001348). The phylogenetic analyses of these sequences determined the circulation of clades 1, 5, and 9 in India. Further, the sequence analysis demonstrated non-synonymous single nucleotide polymorphism (SNPs) among specific ORF of VZV including ORF 14, ORF 22, ORF 36, ORF 37 and ORF 51. Although clade 1 and 5 has been reported earlier, the circulation of clade 9 of VZV has been determined for the first time in India. CONCLUSION: Although the circulation of different clades of VZV was reported from India, the presence of clade 9 was detected for the first time during the mpox surveillance.

Genomic characterization, transcriptome analysis, and pathogenicity of the Nipah virus (Indian isolate).

Nipah virus (NiV) is a high-risk pathogen which can cause fatal infections in humans. The Indian isolate from the 2018 outbreak in the Kerala state of India showed ~ 4% nucleotide and amino acid difference in comparison to the Bangladesh strains of NiV and the substitutions observed were mostly not present in the region of any functional significance except for the phosphoprotein gene. The differential expression of viral genes was observed following infection in Vero (ATCC® CCL-81™) and BHK-21 cells. Intraperitoneal infection in the 10-12-week-old, Syrian hamster model induced dose dependant multisystemic disease characterized by prominent vascular lesions in lungs, brain, kidney and extra vascular lesions in brain and lungs. Congestion, haemorrhages, inflammatory cell infiltration, thrombosis and rarely endothelial syncitial cell formation were seen in the blood vessels. Intranasal infection resulted in respiratory tract infection characterised by pneumonia. The model showed disease characteristics resembling the human NiV infection except that of myocarditis similar to that reported by NiV-Malaysia and NiV-Bangladesh isolates in hamster model. The variation observed in the genome of the Indian isolate at the amino acid levels should be explored further for any functional significance.

Omicron BA.2 lineage predominance in severe acute respiratory syndrome coronavirus 2 positive cases during the third wave in North India.

Background: Recent studies on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reveal that Omicron variant BA.1 and sub-lineages have revived the concern over resistance to antiviral drugs and vaccine-induced immunity. The present study aims to analyze the clinical profile and genome characterization of the SARS-CoV-2 variant in eastern Uttar Pradesh (UP), North India. Methods: Whole-genome sequencing (WGS) was conducted for 146 SARS-CoV-2 samples obtained from individuals who tested coronavirus disease 2019 (COVID-19) positive between the period of 1 January 2022 and 24 February 2022, from three districts of eastern UP. The details regarding clinical and hospitalized status were captured through telephonic interviews after obtaining verbal informed consent. A maximum-likelihood phylogenetic tree was created for evolutionary analysis using MEGA7. Results: The mean age of study participants was 33.9 ± 13.1 years, with 73.5% accounting for male patients. Of the 98 cases contacted by telephone, 30 (30.6%) had a travel history (domestic/international), 16 (16.3%) reported having been infected with COVID-19 in past, 79 (80.6%) had symptoms, and seven had at least one comorbidity. Most of the sequences belonged to the Omicron variant, with BA.1 (6.2%), BA.1.1 (2.7%), BA.1.1.1 (0.7%), BA.1.1.7 (5.5%), BA.1.17.2 (0.7%), BA.1.18 (0.7%), BA.2 (30.8%), BA.2.10 (50.7%), BA.2.12 (0.7%), and B.1.617.2 (1.3%) lineages. BA.1 and BA.1.1 strains possess signature spike mutations S:A67V, S:T95I, S:R346K, S:S371L, S:G446S, S:G496S, S:T547K, S:N856K, and S:L981F, and BA.2 contains S:V213G, S:T376A, and S:D405N. Notably, ins214EPE (S1- N-Terminal domain) mutation was found in a significant number of Omicron BA.1 and sub-lineages. The overall Omicron BA.2 lineage was observed in 79.5% of women and 83.2% of men. Conclusion: The current study showed a predominance of the Omicron BA.2 variant outcompeting the BA.1 over a period in eastern UP. Most of the cases had a breakthrough infection following the recommended two doses of vaccine with four in five cases being symptomatic. There is a need to further explore the immune evasion properties of the Omicron variant.

Genomic profile of SARS-CoV-2 Omicron variant and its correlation with disease severity in Rajasthan.

Background: Omicron, a new variant of Severe Acute Respiratory Syndrome-Coronavirus 2 (SARS-CoV-2), was first detected in November 2021. This was believed to be highly transmissible and was reported to evade immunity. As a result, an urgent need was felt to screen all positive samples so as to rapidly identify Omicron cases and isolate them to prevent the spread of infection. Genomic surveillance of SARS-CoV-2 was planned to correlate disease severity with the genomic profile. Methods: All the SARS-CoV-2 positive cases detected in the state of Rajasthan were sent to our Lab. Samples received from 24 November 2021 to 4 January 2022 were selected for Next-Generation Sequencing (NGS). Processing was done as per protocol on the Ion Torrent S5 System for 1,210 samples and bioinformatics analysis was done. Results: Among the 1,210 samples tested, 762 (62.9%) were Delta/Delta-like and other lineages, 291 (24%) were Omicron, and 157 (12.9%) were invalid or repeat samples. Within a month, the proportion of Delta and other variants was reversed, 6% Omicron became 81%, and Delta and other variants became 19%, initially all Omicron cases were seen in international travelers and their contacts but soon community transmission was seen. The majority of patients with Omicron were asymptomatic (56.7%) or had mild disease (33%), 9.2% had moderate symptoms, and two (0.7%) had severe disease requiring hospitalization, of which one (0.3%) died and the rest were (99.7%) recovered. History of vaccination was seen in 81.1%, of the previous infection in 43.2% of cases. Among the Omicron cases, BA.1 (62.8%) was the predominant lineage followed by BA.2 (23.7%) and B.1.529 (13.4%), rising trends were seen initially for BA.1 and later for BA.2 also. Although 8.9% of patients with Delta lineage during that period were hospitalized, 7.2% required oxygen, and 0.9% died. To conclude, the community spread of Omicron occurred in a short time and became the predominant circulating variant; BA.1 was the predominant lineage detected. Most of the cases with Omicron were asymptomatic or had mild disease, and the mortality rate was very low as compared to Delta and other lineages.

Delta variant SARS-CoV-2 infections in pediatric cases during the second wave in India.

BACKGROUND: During October 2020, Delta variant was detected for the first time in India and rampantly spread across the globe. It also led to second wave of pandemic in India which affected millions of people. However, there is limited information pertaining to the SARS-CoV-2 strain infecting the children in India. METHODS: Here, we assessed the SARS-CoV-2 lineages circulating in the pediatric population of India during the second wave of the pandemic. Clinical and demographic details linked with the nasopharyngeal/oropharyngeal swabs (NPS/OPS) collected from SARS-CoV-2 cases (n = 583) aged 0-18 year and tested positive by real-time RT-PCR were retrieved from March to June 2021. RESULTS: Symptoms were reported among 37.2% of patients and 14.8% reported to be hospitalized. The E gene CT value had significant statistical difference at the point of sample collection when compared to that observed in the sequencing laboratory. Out of these 512 sequences 372 were VOCs, 51 were VOIs. Most common lineages observed were Delta, followed by Kappa, Alpha and B.1.36, seen in 65.82%, 9.96%, 6.83% and 4.68%, respectively in the study population. CONCLUSION: Overall, it was observed that Delta strain was the leading cause of SARS-CoV-2 infection in Indian children during the second wave of the pandemic. We emphasize on the need of continuous genomic surveillance in SARS-CoV-2 infection even amongst children.

Zika a Vector Borne Disease Detected in Newer States of India Amidst the COVID-19 Pandemic.

Background: During the second wave of the COVID-19 pandemic, outbreaks of Zika were reported from Kerala, Uttar Pradesh, and Maharashtra, India in 2021. The Dengue and Chikungunya negative samples were retrospectively screened to determine the presence of the Zika virus from different geographical regions of India. Methods: During May to October 2021, the clinical samples of 1475 patients, across 13 states and a union territory of India were screened and re-tested for Dengue, Chikungunya and Zika by CDC Trioplex Real time RT-PCR. The Zika rRTPCR positive samples were further screened with anti-Zika IgM and Plaque Reduction Neutralization Test. Next generation sequencing was used for further molecular characterization. Results: The positivity was observed for Zika (67), Dengue (121), and Chikungunya (10) amongst screened cases. The co-infections of Dengue/Chikungunya, Dengue/Zika, and Dengue/Chikungunya/Zika were also observed. All Zika cases were symptomatic with fever (84%) and rash (78%) as major presenting symptoms. Of them, four patients had respiratory distress, one presented with seizures, and one with suspected microcephaly at birth. The Asian Lineage of Zika and all four serotypes of Dengue were found in circulation. Conclusion: Our study indicates the spread of the Zika virus to several states of India and an urgent need to strengthen its surveillance.

First case of Zika virus infection during an outbreak of chikungunya in a rural region of Maharashtra state, India.

BACKGROUND: In July 2021, an outbreak of chikungunya virus (CHIKV) was reported in a rural region of Maharashtra state, India. METHODS: Serum samples of symptomatic cases (n=33) were screened for dengue virus (DENV), CHIKV and Zika virus (ZIKV) by molecular and serological assays. RESULTS: The first case of ZIKV infection from Maharashtra was detected and confirmed by molecular and serological assays. Complete genome sequencing revealed that the ZIKV sequence belongs to the Asian genotype and had a closer homology with pre-epidemic strains present before 2007. CONCLUSIONS: ZIKV surveillance needs to be strengthened in the regions experiencing dengue and chikungunya outbreaks.

Nipah Virus Outbreak in Kerala State, India Amidst of COVID-19 Pandemic.

We report here a Nipah virus (NiV) outbreak in Kozhikode district of Kerala state, India, which had caused fatal encephalitis in a 12-year-old boy and the outbreak response, which led to the successful containment of the disease and the related investigations. Quantitative real-time reverse transcription (RT)-PCR, ELISA-based antibody detection, and whole genome sequencing (WGS) were performed to confirm the NiV infection. Contacts of the index case were traced and isolated based on risk categorization. Bats from the areas near the epicenter of the outbreak were sampled for throat swabs, rectal swabs, and blood samples for NiV screening by real-time RT-PCR and anti-NiV bat immunoglobulin G (IgG) ELISA. A plaque reduction neutralization test was performed for the detection of neutralizing antibodies. Nipah viral RNA could be detected from blood, bronchial wash, endotracheal (ET) secretion, and cerebrospinal fluid (CSF) and anti-NiV immunoglobulin M (IgM) antibodies from the serum sample of the index case. Rapid establishment of an onsite NiV diagnostic facility and contact tracing helped in quick containment of the outbreak. NiV sequences retrieved from the clinical specimen of the index case formed a sub-cluster with the earlier reported Nipah I genotype sequences from India with more than 95% similarity. Anti-NiV IgG positivity could be detected in 21% of Pteropus medius (P. medius) and 37.73% of Rousettus leschenaultia (R. leschenaultia). Neutralizing antibodies against NiV could be detected in P. medius. Stringent surveillance and awareness campaigns need to be implemented in the area to reduce human-bat interactions and minimize spillover events, which can lead to sporadic outbreaks of NiV.

Detection and isolation of SARS-CoV-2 Eta variant from the international travelers and local residents of India.

International travel has been the major source for the rapid spread of new SARS-CoV-2 variants across the globe. During SARS-CoV-2 genomic surveillance, a total of 212 SARS-CoV-2 positive clinical specimens were sequenced using next-generation sequencing. A complete SARS-CoV-2 genome could be retrieved from 90 clinical specimens. Of them, 14 sequences belonged to the Eta variant from clinical specimens of international travelers (n = 12) and local residents (n = 2) of India, and 76 belonged to other SARS-CoV-2 variants. Of all the Eta-positive specimens, the virus isolates were obtained from the clinical specimens of six international travelers. Many variants of interest have been found to cause substantial community transmission or cluster infections. The detection of this variant with lethal E484K mutation across the globe and India necessitates persistent genomic surveillance of the SARS-CoV-2 variants, which would aid in taking preventive action.

Isolation and characterization of SARS-CoV-2 Beta variant from UAE travelers.

BACKGROUND: The emergence of SARS-CoV-2 variants in places where the virus is uncontained poses a global threat from the perspective of public health and vaccine efficacy. Travel has been important factor for the easy spread of SARS-CoV-2 variants worldwide. India has also observed the importation of SARS-CoV-2 variants through international travelers. METHODS: In this study, we have collected the oropharyngeal and nasopharyngeal swab specimens from 58 individuals with travel history from United Arab Emirates (UAE), East, West and South Africa, Qatar, Ukraine and Saudi Arabia arrived in India during February-March 2021. The clinical specimens were initially screened for SARS-CoV-2 using Real time RT-PCR. All the specimens were inoculated on to Vero CCL-81 cells for virus isolation. The viral isolates were further sequenced using Next-Generation Sequencing. RESULTS: All 58 cases were tested positive for SARS-CoV-2 using Real time RT-PCR. Four specimens showed progressive infectivity with fusion of the infected cells with neighboring cells leading to large mass of cells. Replication competent virus was confirmed from culture supernatant of the passage 2 using Real time RT-PCR. Two plaque purified SARS-CoV-2 isolates demonstrated high viral RNA load of 3.8-7.5 × 1011 and 1.1-1.6 × 1011 at passage 4 and 5 respectively. Nucleotide variations along with amino acid changes were also observed among these two isolates at passage 2-5. All four cases were male with no symptoms and co-morbidity. The sequence analysis has shown two different clusters, first cluster with nucleotide deletions in the ORF1ab and the spike, while second cluster with deletions in spike region. The viral isolates demonstrated 99.88-99.96% nucleotide identity with the representative sequences of Beta variant (B.1.351). CONCLUSION: These findings suggest easier transmission of SARS-CoV-2 variants with human mobility through international travel. The isolated Beta variant would be useful to determine the protective efficacy of the currently available and upcoming COVID-19 vaccines in India.

Possible Role of Accessory Proteins in the Viral Replication for the 20I/501Y.V1 (B.1.1.7) SARS CoV-2 Variant.

The emergence of new severe acute respiratory syndrome coronavirus-2 (SARS CoV-2) has been a global concern. The B.1.1.7 variant of SARS CoV-2 is reported to cause higher transmission. The study investigates the replication cycle and transcriptional pattern of the B.1.1.7 to hypothesis the possible role of different genes in viral replication. It was observed that the B.1.1.7 variant required a longer maturation time. The transcriptional response demonstrated higher expression of ORF6 and ORF8 compared to nucleocapsid transcript till the eclipse period which might influence higher viral replication. The number of infectious viruses titer is higher in the B.1.1.7, despite a lesser copy number than B.1, indicating higher transmissibility. The experimental evidence published linked ORF6 and ORF8 to play important role in replication and we also observed their higher expression. This leads us to hypothesis the possible role of ORF6 and ORF8 in B.1.1.7 higher replication which causes higher transmission.

Clinical Characterization and Genomic Analysis of Samples from COVID-19 Breakthrough Infections during the Second Wave among the Various States of India.

From March to June 2021, India experienced a deadly second wave of COVID-19, with an increased number of post-vaccination breakthrough infections reported across the country. To understand the possible reason for these breakthroughs, we collected 677 clinical samples (throat swab/nasal swabs) of individuals from 17 states/Union Territories of the country who had received two doses (n = 592) and one dose (n = 85) of vaccines and tested positive for COVID-19. These cases were telephonically interviewed and clinical data were analyzed. A total of 511 SARS-CoV-2 genomes were recovered with genome coverage of higher than 98% from both groups. Analysis of both groups determined that 86.69% (n = 443) of them belonged to the Delta variant, along with Alpha, Kappa, Delta AY.1, and Delta AY.2. The Delta variant clustered into four distinct sub-lineages. Sub-lineage I had mutations in ORF1ab A1306S, P2046L, P2287S, V2930L, T3255I, T3446A, G5063S, P5401L, and A6319V, and in N G215C; Sub-lineage II had mutations in ORF1ab P309L, A3209V, V3718A, G5063S, P5401L, and ORF7a L116F; Sub-lineage III had mutations in ORF1ab A3209V, V3718A, T3750I, G5063S, and P5401L and in spike A222V; Sub-lineage IV had mutations in ORF1ab P309L, D2980N, and F3138S and spike K77T. This study indicates that majority of the breakthrough COVID-19 clinical cases were infected with the Delta variant, and only 9.8% cases required hospitalization, while fatality was observed in only 0.4% cases. This clearly suggests that the vaccination does provide reduction in hospital admission and mortality.

An Epidemiological Analysis of SARS-CoV-2 Genomic Sequences from Different Regions of India.

The number of Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) cases is increasing in India. This study looks upon the geographic distribution of the virus clades and variants circulating in different parts of India between January and August 2020. The NPS/OPS from representative positive cases from different states and union territories in India were collected every month through the VRDLs in the country and analyzed using next-generation sequencing. Epidemiological analysis of the 689 SARS-CoV-2 clinical samples revealed GH and GR to be the predominant clades circulating in different states in India. The northern part of India largely reported the 'GH' clade, whereas the southern part reported the 'GR', with a few exceptions. These sequences also revealed the presence of single independent mutations-E484Q and N440K-from Maharashtra (first observed in March 2020) and Southern Indian States (first observed in May 2020), respectively. Furthermore, this study indicates that the SARS-CoV-2 variant (VOC, VUI, variant of high consequence and double mutant) was not observed during the early phase of virus transmission (January-August). This increased number of variations observed within a short timeframe across the globe suggests virus evolution, which can be a step towards enhanced host adaptation.

Immunogenicity and protective efficacy of inactivated SARS-CoV-2 vaccine candidate, BBV152 in rhesus macaques.

The COVID-19 pandemic is a global health crisis that poses a great challenge to the public health system of affected countries. Safe and effective vaccines are needed to overcome this crisis. Here, we develop and assess the protective efficacy and immunogenicity of an inactivated SARS-CoV-2 vaccine in rhesus macaques. Twenty macaques were divided into four groups of five animals each. One group was administered a placebo, while three groups were immunized with three different vaccine candidates of BBV152 at 0 and 14 days. All the macaques were challenged with SARS-CoV-2 fourteen days after the second dose. The protective response was observed with increasing SARS-CoV-2 specific IgG and neutralizing antibody titers from 3rd-week post-immunization. Viral clearance was observed from bronchoalveolar lavage fluid, nasal swab, throat swab and lung tissues at 7 days post-infection in the vaccinated groups. No evidence of pneumonia was observed by histopathological examination in vaccinated groups, unlike the placebo group which exhibited interstitial pneumonia and localization of viral antigen in the alveolar epithelium and macrophages by immunohistochemistry. This vaccine candidate BBV152 has completed Phase I/II (NCT04471519) clinical trials in India and is presently in phase III, data of this study substantiates the immunogenicity and protective efficacy of the vaccine candidates.

Genome Sequencing Reveals a Mixed Picture of SARS-CoV-2 Variant of Concern Circulation in Eastern Uttar Pradesh, India.

Uttar Pradesh is the densely populated state of India and is the sixth highest COVID-19 affected state with 22,904 deaths recorded on November 12, 2021. Whole-genome sequencing (WGS) is being used as a potential approach to investigate genomic evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. In this study, a total of 87 SARS-CoV-2 genomes-49 genomes from the first wave (March 2020 to February 2021) and 38 genomes from the second wave (March 2021 to July 2021) from Eastern Uttar Pradesh (E-UP) were sequenced and analyzed to understand its evolutionary pattern and variants against publicaly available sequences. The complete genome analysis of SARS-CoV-2 during the first wave in E-UP largely reported transmission of G, GR, and GH clades with specific mutations. In contrast, variants of concerns (VOCs) such as Delta (71.0%) followed by Delta AY.1 (21.05%) and Kappa (7.9%) lineages belong to G clade with prominent signature amino acids were introduced in the second wave. Signature substitution at positions S:L452R, S:P681R, and S:D614G were commonly detected in the Delta, Delta AY.1, and Kappa variants whereas S:T19R and S:T478K were confined to Delta and Delta AY.1 variants only. Vaccine breakthrough infections showed unique mutational changes at position S:D574Y in the case of the Delta variant, whereas position S:T95 was conserved among Kappa variants compared to the Wuhan isolate. During the transition from the first to second waves, a shift in the predominant clade from GH to G clade was observed. The identified spike protein mutations in the SARS-CoV-2 genome could be used as the potential target for vaccine and drug development to combat the effects of the COVID-19 disease.